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Abstract 

A first measurement is presented of the charge asymmetry in the hadronic final state 
from the hard interaction in deep-inelastic ep neutral current scattering at HERA. The mea- 
surement is performed in the range of negative squared four momentum transfer 100 < 
Q 2 < 8,000 GeV 2 . The difference between the event normalised distributions of the 
scaled momentum, x p , for positively and negatively charged particles, measured in the cur- 
rent region of the Breit frame, is studied together with its evolution as a function of Q. The 
results are compared to Monte Carlo models at the hadron and parton level. 
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1 Introduction 



In lepton proton deep-inelastic scattering (DIS) at large Bjorken x the contribution of u valence 
quarks from the proton to the hard interaction dominates over that from the d valence quarks due 
to their larger charge and greater abundance. Hence an asymmetry in the number of positively 
and negatively charged particles is observed in the final state [1]. It has been demonstrated that 
the charge sign asymmetry of the hadronic final state in pp collisions at RHIC [2] is sensitive to 
the valence quark distribution [3]. 

In a recent paper, HI presented a study of the inclusive charged particle production in high 
Q 2 deep-inelastic scattering at HERA [4]. The measurement is performed in the current hemi- 
sphere of the Breit frame [5]. In the naive quark parton model (QPM) the momentum of the 
scattered parton in the Breit frame is Q/2, where Q 2 is the virtuality of the exchanged boson. 
The main observable is x p , the charged particle momentum in the current region of the Breit 
frame scaled to Q/2. General agreement was observed between ep, e + e~ data and Monte Carlo 
predictions, broadly supporting the concept of quark fragmentation universality. Hadrons with 
small values of x p are predominately produced by fragmentation, while hadrons at large x p are 
more likely to contain a parton from the hard interaction. Therefore a study of the x v distribu- 
tion separately for positively and negatively charged particles should reveal information about 
the valence quarks and their fragmentation. 

The analysis presented here utilises the same data and methodology as in [4] but separates 
the positively and negatively charged particles into different distributions. In addition the charge 
asymmetry is studied. The results are compared with predictions from different fragmentation 
models implemented in Monte Carlo programs. 



2 Experimental Method 

A full description of the HI detector can be found elsewhere [6] and only those components 
most relevant for this analysis are mentioned briefly here. The origin of the HI coordinate 
system is the nominal ep interaction point, the direction of the proton beam defining the positive 
z-axis (forward region). 

The Liquid Argon (LAr) calorimeter measures the positions and energies of particles, in- 
cluding that of the scattered positron, over the polar angle range 4° < 9 < 154°. The calorime- 
ter consists of an electromagnetic section with lead absorbers and a hadronic section with steel 
absorbers. The energy resolution for electrons in the electromagnetic section is a{E)/E = 
11.5%/VE [GeV] © 1% [7]. 

Charged particles are measured in the Central Tracking Detector (CTD) in the range 
20° < 9 < 165°. The CTD comprises two large cylindrical Central Jet Chambers (CJCs) 
arranged concentrically around the beam-line, complemented by a silicon vertex detector [8] 
covering the range 30° < 9 < 150°, two z-drift chambers and two multiwire proportional 
chambers for triggering purposes, all within a solenoidal magnetic field of strength 1.16 T. The 
transverse momentum resolution is <j(pt) /pt — 0.006 pt [GeV] © 0.02 [9] . In each event the 
tracks are used in a common fit procedure to determine the ep interaction vertex. 
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The data used in this analysis correspond to an integrated luminosity of 44 pb" 1 and were 
taken by HI in the year 2000 when protons with an energy of 920 GeV collided with positrons 
with an energy of 27.5 GeV. The event and track selection follows that of [4], here only the 
kinematic phase space is defined. 

The scattered positron is detected in the LAr calorimeter in the polar angular range 10° < 
9 e < 150° and with energy greater than 11 GeV. The negative squared four momentum transfer 
is required to be in the range 100 < Q 2 < 8 000 GeV 2 and the inelasticity y, which is the 
fractional energy loss of the positron in the proton rest frame, to be in the range 0.05 < y < 0.6. 
The polar scattering angle for a massless parton, calculated from the positron kinematics in the 
quark-parton model (QPM) approximatiorQ, is required to be in the range 30° < 9 q i ab < 150°. 
This ensures that the current region of the Breit frame remains in the central region of the 
detector where there is high acceptance and track reconstruction efficiency. It should be noted 
that the kinematic phase space is defined solely from the scattered electron and can be applied 
in a simple way to theoretical models. The full event selection outlined above results in a data 
sample of about 60,000 events. 

3 Observables 

The current hemisphere of the Breit frame of reference provides a kinematic region where the 
properties of the scattered quark can be studied with a well defined and relatively clean separa- 
tion from the proton remnants. In the Breit frame the virtual space-like photon has momentum 
Q but no energy. The photon direction defines the negative z'-axis and the current hemisphere. 
Within the QPM the photon collides head on with a massless quark of longitudinal momentum 
Q/2. The struck quark thus scatters with an equal but opposite momentum into the current 
hemisphere while the proton remnants go into the opposite (target) hemisphere. The energy 
scale is set by the virtual photon at Q/2. The boost to the Breit frame is defined using kine- 
matics calculated from the scattered positron and is thus not biased by mismeasurements of the 
hadronic final state. Hadrons emerging from the interaction with negative longitudinal momenta 
in this frame are assigned to the current region and associated with the struck quark. 

The scaled momentum variable, x p , is defined to be pf l /(Q/2) where ph is the momentum 
of a charged track in the current region of the Breit frame. The inclusive, event normalised, 
charged particle scaled momentum distribution, D(x p , Q), is calculated as Tf-pS where in each 
Q range, N is the number of selected events and dn is the total number of charged tracks 
with scaled momentum x p in the interval dx p . The distribution D(x p , Q) has been measured 
previously [4]. 

In the present study, the scaled momentum distribution is defined separately for positively 
charged particles, D + (x p ,Q), and for negatively charged particles, D~(x p ,Q). The charge 
asymmetry is defined as A(x p , Q) = (D + (x p , Q) — D~(x p , Q))/D(x p , Q), where the relation 
D(x p ,Q) = D + (x p ,Q) + D~(x p ,Q) holds. Systematic errors mostly cancel in the charge 
asymmetry. 

'in the approximation of a massless scattered quark, the polar scattering angle is given by 
0q,iab = cos -1 ((xs(xs — Q 2 ) — AE 2 Q 2 ) / (xs(xs — Q 2 ) + AE 2 Q 2 )), where E is the incoming positron beam 
energy, s is the ep centre of mass energy squared and Bjorken x is the fraction of the proton momentum carried by 
the struck quark in the QPM. 
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4 Phenomenology 



A more complete discussion on the different models of the parton cascade and hadronisation 
processes can be found in [4]. 

The data are compared to predictions of the Parton Shower model (PS) [10], as imple- 
mented in the RAPGAP [11] Monte Carlo program and to predictions of the Colour Dipole 
Model (CDM) [12] both matched to 0(a s ) matrix elements. ARIADNE [13] provides an im- 
plementation of CDM and is used in the DJANGO [14] Monte Carlo program. Both the PS 
and CDM predictions use the Lund string model for hadronisation [15]. The HERWIG Monte 
Carlo [16] program uses the parton shower model to describe the fragmentation process but 
incorporates the cluster model of hadronisation [17]. 

In the Soft Colour Interaction model (SCI) [18] soft gluons are exchanged between the 
partons produced in the parton shower and the proton remnant. A refined version of the model 
uses a generalised area law (GAL) [19] for the colour rearrangement probability. Predictions 
for SCI and GAL models are obtained using the LO generator programs LEPTO [20] for DIS. 
Higher order QCD effects are simulated using parton showers. Hadronisation is simulated using 
the Lund string model. The SCI and GAL models produce similar predictions for the charge 
asymmetry, so only GAL model predictions are given in this paper. 

It is possible to turn off the hadronisation and compare data with parton level predictions 
using the assumption of local parton hadron duality. This has been done with the CDM predic- 
tions where a quark, with fractional charge, is taken as equivalent to a charged hadron of unit 
charge. The predictions are made after the main parton cascade has taken place and the gluons 
are ignored. 

The CTEQ5L [21] parton density function (PDF) is used for all model predictions. Other 
PDFs [22, 23] lead to charge asymmetries in agreement with the prediction based on CTEQ5L 
to within ± 0.01. 



5 Data Correction 

The data are corrected for detector acceptance, efficiency and resolution effects using Monte 
Carlo event samples generated with the RAPGAP and DJANGO programs. All generated events 
are passed through the GEANT [24] based simulation of the HI apparatus and are reconstructed 
and analysed using the same programs as used for the data. These Monte Carlo event samples 
give a good description of the data. The residual contribution of charged particles from the 
weak decay of neutral particles (e.g. K° and A's) is about 8% and is subtracted from the data as 
part of the correction procedure. The effects of QED radiation are corrected for using the HER- 
ACLES [14] program incorporated within the above Monte Carlos. The total correction factor 
is calculated using DJANGO from the ratio of the number of entries in each bin at hadron level 
to that at detector level. The bin sizes are chosen to give high acceptance and purit^H, typically 

2 The acceptance (purity) is defined as the ratio of the number of charged hadrons generated and reconstructed 
to the total number of charged hadrons generated (reconstructed) in that bin. 
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above 60% with a minimum of 40%. The total correction factor applied to the uncorrected data 
points is typically 1.0 — 1.2 for D ± (x p , Q). In general the uncertainty in the boost to the Breit 
frame dominates the resolution in x p . The correction associated with the tracking dominates in 
the highest Q 2 region where there is a somewhat reduced acceptance for the current region of 
the Breit frame within the CTD. 

In the measurement of the charge asymmetry A(x p , Q) contributions to the correction factor 
such as QED corrections, efficiencies and acceptance mostly cancel and as a result the correc- 
tion factor is consistent with 1.0. 

The data correction method was cross checked using a matrix migration unfolding proce- 
dure, which was found to be in agreement. 



6 Systematic Uncertainties 

The positron energy scale uncertainty is 0.7 — 3 % depending on the position of the detected 
positron in the LAr calorimeter. This uncertainty affects both the phase space and boost calcula- 
tion. The resulting uncertainty on D ± (x p , Q) is independent of Q but varies with x p from 0.5% 
(x p ~ 0.1) to 11% (x p ~ 1.0). The precision in the reconstruction of the scattered positron 
direction leads to a systematic error of about 1%. The hadronic energy scale uncertainty is 
4% and this leads to an error of less than 1% on D ± (x p , Q). The uncertainty in the correction 
factor arising from using different Monte Carlo models in the correction procedure, taken as 
the full difference between correcting the data with RAPGAP or DJANGO, results in a typical 
error of 1.5% on D ± (x p , Q). In the defined kinematic region, errors arising from non ep back- 
ground are negligible. The systematic error associated with the track reconstruction (e.g. track 
reconstruction efficiency, vertex reconstruction efficiency, weak decays and nuclear interaction 
uncertainties) is estimated to be 2.5% for D ± (x p , Q). All sources of error are treated as un- 
correlated, apart from the positron energy scale uncertainty which is treated as fully correlated 
between bins. 

For the charge asymmetry measurements the effects of the systematic errors mostly cancel 
and the only significant contribution is from the track reconstruction uncertainty. The track 
reconstruction efficiency for positive and negative tracks is conservatively considered as fully 
anti correlated in the estimation of the effect on the charge asymmetry. Secondary interactions 
of tracks in the material in front of the trackers may also lead to charge asymmetries. At large 
track momentum a possible charge bias in the track reconstruction is studied using the track 
of the scattered lepton in e + p and e~p data (in the year 1998-1999 HERA operated with an 
electron beam). These sources of systematic error result in an uncertainty of 2.5% on the charge 
asymmetry. 



7 Results 

The scaled momentum distribution for all charged particles, and for positively, and negatively 
charged particles separately, is shown in figure \T^l) and table [Q There are significantly more 
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particles produced at low x p than at high x p , in agreement with predictions. The scaled momen- 
tum distribution for positively and negatively charged particles are very similar at low x p but at 
high x p there is a clear excess of positively charged particles. 

The charge asymmetry can be as large as 0.18 as shown in figure dh) and table [21 The scaled 
momentum distribution and its asymmetry is described by the Monte Carlo models. Models 
incorporating string hadronisation (PS, CDM, GAL) produce a smaller charge asymmetry at 
high x p than that produced by the cluster hadronisation model (HERWIG). Similar differences 
between the models at high x p were observed in [4]. 

In figure dh) the charge asymmetry is compared to predictions from CDM before and after 
hadronisation. It is observed that at high x p the hadron and quark levels are in good agreement 
and both agree with the data. However, as x p gets smaller a large difference develops with the 
quark level asymmetry prediction constant at about 0.12, while the hadron level, and the data, 
fall to zero. This is consistent with the expectation that the hadrons at low x p are dominantly 
produced by fragmentation while hadrons at high x p retain the memory of the charge of the 
scattered quark from the hard interaction. It should be noted that sea quarks and gluons will 
produce, on average, charge symmetric hadronic final states, reducing the charge asymmetry 
expected from valence quarks alone. 

The analysis intervals in Q 2 and the average values of Q and Bjorken x are shown in table[3] 
Figures [2] and |3] (tables H] and [5]) show the scaled momentum distribution and charge asymmetry 
as a function of Q in different x p intervals. The charge asymmetry observed at large x p evolves 
to larger values as Q increases. The largest asymmetries of about 0.4 are obtained in the highest 
Q and highest x p intervals. It should be noted that higher average Q corresponds to higher 
average Bjorken x (table [3]) and hence the highest Q intervals are most sensitive to the valence 
quark distribution. The scaled momentum distributions (D ± (x p , Q)) are broadly predicted by 
the Monte Carlo predictions but tend to undershoot at large Q, similar to their sum [4]. The 
charge asymmetry is well described by the Monte Carlo. 

8 Conclusions 

The first measurement of the charge asymmetry of the hadronic final state at HERA is pre- 
sented. The charge asymmetry is found to be dependent on the scaled momentum x p with a 
larger asymmetry for large x p . The observed charge asymmetry at large x p is found to increase 
with the scale Q corresponding at HERA to an enhancement at large Bjorken x. The results 
are consistent with the expectation that at high x p the asymmetry is directly related to the va- 
lence quark content of the proton. The observed charge asymmetry is reproduced by various 
models. The data are expected to provide useful information for the extraction of fragmentation 
functions and additional constraints on the valence quark distribution of the proton. 
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Xp 


£> (±) (* P ) 


8stat [%. 


Stot[%] 


8scale[%_ 




sum 


0.0 < x p < 0.02 


17.26 


0.8 


2.8 


0.7/0.5 


0.02 <x„< 0.05 


33.55 


0.5 


2.6 


1.1/1.2 


0.05 <x p < 0.1 


23.35 


0.4 


2.6 


1.3/1.3 


0.1 < x p < 0.2 


10.27 


0.5 


2.6 


1.4/1.6 


0.2 < x p < 0.3 


3.98 


0.7 


2.9 


2.0/1.9 


0.3 < x p < 0.4 


1.766 


1.1 


3.2 


2.6/2.9 


0.4 < x p < 0.5 


0.843 


1.6 


3.6 


3.8/3.4 


0.5 < x p < 0.7 


0.326 


1.8 


3.6 


5.0/5.1 


0.7 <x v < 1.0 


0.0553 


3.4 


5.5 


9.0/8.7 




positive 


0.0 < a; p < 0.02 


8.74 


1.1 


3.0 


0.7/0.4 


0.02 < Xp < 0.05 


16.87 


0.7 


2.7 


1.1/1.1 


0.05 < Xp < 0.1 


11.70 


0.7 


2.7 


1.3/1.3 


0.1 < Xp < 0.2 


5.18 


0.7 


2.7 


1.3/1.4 


0.2 < Xp < 0.3 


2.05 


1.0 


2.9 


1.8/1.6 


0.3 < x p < 0.4 


0.930 


1.5 


3.4 


2.7/3.1 


0.4 < x p < 0.5 


0.460 


2.2 


3.9 


3.9/3.3 


0.5 < Xp < 0.7 


0.180 


2.5 


4.0 


4.1/5.0 


0.7 < x p < 1.0 


0.0327 


4.6 


6.3 


9.0/9.0 




negative 


0.0 < x p < 0.02 


8.52 


1.2 


3.0 


0.7/0.5 


0.02 < Xp < 0.05 


16.68 


0.7 


2.7 


1.2/1.3 


0.05 < x p < 0.1 


11.66 


0.6 


2.7 


1.3/1.3 


0.1 < x p < 0.2 


5.09 


0.7 


2.7 


1.4/1.7 


0.2 < x v < 0.3 


1.93 


1.1 


2.9 


2.2/2.2 


0.3 < Xp < 0.4 


0.836 


1.6 


3.4 


2.6/2.7 


0.4 < x p < 0.5 


0.383 


2.4 


4.0 


3.7/3.4 


0.5 < x p < 0.7 


0.146 


2.7 


4.1 


5.6/5.1 


0.7 < x p < 1.0 


0.0226 


5.4 


6.9 


9.0/11.0 



Table 1 : The measured normalised distribution of the scaled momentum for all charged parti- 
cles, D(x p ), and for positively, D + (x p ), and negatively, D~(x p ), charged particles, for different 
x p intervals, shown with the statistical error (5 sta t), the total error including statistical and sys- 
tematic errors added in quadrature (S to t ), and the correlated error coming from the electron 
energy scale uncertainty (S sca i e ) which is shown as two errors (+/— ) and not included in the 
total error. The average Q value for the data is 19.5 GeV. 
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Xp 


A(x p ) 


$stat 


Stat 


0.0 < x p < 0.02 


0.013 


0.008 


0.027 


0.02 < x p < 0.05 


0.006 


0.005 


0.026 


0.05 < Xp < 0.1 


0.002 


0.004 


0.026 


0.1 < x p < 0.2 


0.009 


0.005 


0.026 


0.2 < Xp < 0.3 


0.030 


0.007 


0.027 


0.3 < x p < 0.4 


0.054 


0.011 


0.029 


0.4 < x p < 0.5 


0.087 


0.016 


0.033 


0.5 < x p < 0.7 


0.106 


0.018 


0.035 


0.7 < x p < 1.0 


0.181 


0.034 


0.047 



Table 2: The charge asymmetry, A(x p ), for different x p intervals, shown with the statistical 
error (5 stat ), the total error including statistical and systematic errors added in quadrature (S tot ). 
The average Q value for the data is 19.5 GeV. 



Q 2 [GeV 2 ] 


(Q) [GeV] 


5Q [GeV] 


(x) 


5x 


100 < Q 2 < 175 


12.3 


0.1 


0.00370 


0.00004 


175 < Q 2 < 250 


14.5 


0.1 


0.00952 


0.00007 


250 < Q 2 < 450 


18.0 


0.1 


0.1559 


0.0001 


450 < Q 2 < 1000 


25.0 


0.3 


0.0254 


0.0003 


1000 < Q 2 < 2000 


36.6 


0.8 


0.044 


0.001 


2000 < Q 2 < 8000 


58.5 


2.1 


0.087 


0.003 



Table 3: Average Q and x values and their statistical errors for the selected events in the Q 2 
intervals used in this analysis. The average Q value for all data is 19.5 GeV. 
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positive 


negative 


Q 2 [GcY 2 ] 


D+(x p ,Q) | 6 stat [%] | 5 4o t[%] | ■5 8eo!e [%] 


D-(x p ,Q) | S stat l*\ | Stat [%] | S scale l%] 




0.0 < x p < 0.02 


100 < Q 2 < 175 


3.14 


5.2 


6.5 


1.6/0.3 


3.08 


5.3 


6.6 


1.5/0.8 


175 < Q 2 < 250 


4.79 


2.7 


3.9 


0.8/0.6 


4.75 


2.7 


4.0 


0.7/0.7 


250 < Q 2 < 450 


7.87 


2.1 


3.5 


1.0/0.3 


7.71 


2.1 


3.4 


0.3/0.3 


450 < Q 2 < 1000 


15.39 


2.3 


3.6 


0.8/0.4 


14.56 


2.3 


3.5 


0.3/0.8 


1000 < Q 2 < 2000 


28.00 


3.2 


4.4 


0.2/1.0 


28.62 


3.2 


4.1 


0.4/0.8 


2000 < Q 2 < 8000 


53.58 


3.8 


4.6 


0.7/0.7 


50.23 


3.9 


4.8 


0.8/1.3 




0.02 < a 


p < 0.05 


100 < Q 2 < 175 


9.52 


2.2 


4.0 


1.4/0.8 


9.54 


2.2 


3.4 


1.2/1.4 


175 < Q 2 < 250 


13.29 


1.3 


2.9 


1.0/1.2 


12.97 


1.3 


2.9 


1.2/1.7 


250 < Q 2 < 450 


17.52 


1.1 


2.9 


1.0/1.2 


17.72 


1.1 


3.0 


1.0/1.1 


450 < Q 2 < 1000 


24.84 


1.5 


3.0 


0.8/1.4 


24.14 


1.5 


3.0 


1.2/1.4 


1000 < Q 2 < 2000 


31.55 


2.5 


3.6 


0.5/1.0 


30.64 


2.5 


4.5 


0.8/1.1 


2000 < Q 2 < 8000 


38.96 


3.8 


5.9 


2.2/0.3 


37.02 


3.9 


5.1 


1.2/1.0 




0.05 < 


Up < 0.1 


100 < Q 2 < 175 


9.01 


1.8 


3.2 


1.6/1.1 


8.81 


1.8 


3.3 


1.5/0.8 


175 < Q 2 < 250 


10.89 


1.1 


2.9 


1.1/1.6 


10.91 


1.1 


3.0 


1.0/1.3 


250 < Q 2 < 450 


12.46 


1.0 


2.8 


1.3/1.2 


12.45 


1.0 


2.8 


1.4/1.4 


450 < Q 2 < 1000 


13.55 


1.5 


3.0 


1.7/1.2 


13.47 


1.6 


3.2 


1.4/1.1 


1000 < Q 2 < 2000 


14.79 


2.8 


3.6 


0.7/1.1 


14.90 


2.9 


3.9 


0.6/1.4 


2000 < Q 2 < 8000 


13.91 


4.9 


5.5 


1.2/1.3 


13.37 


5.0 


7.4 


1.3/1.5 




0.1 < a 


p < 0.2 


100 < Q 2 < 175 


4.69 


1.7 


3.2 


1.1/1.1 


4.68 


1.7 


3.3 


1.2/0.8 


175 < Q 2 < 250 


5.06 


1.1 


2.8 


1.7/1.5 


5.08 


1.1 


2.9 


1.6/2.0 


250 < Q 2 < 450 


5.37 


1.1 


3.0 


1.2/1.6 


5.17 


1.1 


3.1 


1.5/2.0 


450 < Q 2 < 1000 


5.44 


1.7 


3.0 


1.2/1.0 


5.39 


1.8 


3.1 


1.2/1.4 


1000 < Q 2 < 2000 


5.55 


3.3 


3.8 


1.0/1.4 


5.08 


3.5 


6.0 


1.0/1.2 


2000 < Q 2 < 8000 


5.16 


5.8 


6.7 


1.6/2.3 


4.67 


6.2 


7.3 


1.1/2.3 




0.2 < a 


p < 0.3 


100 < Q 2 < 175 


1.95 


2.7 


4.1 


1.6/1.2 


1.86 


2.8 


4.0 


2.0/1.8 


175 < Q 2 < 250 


2.08 


1.8 


3.4 


2.0/1.8 


1.96 


1.8 


3.4 


2.7/2.1 


250 < Q 2 < 450 


2.09 


1.8 


3.3 


1.7/2.2 


1.98 


1.8 


3.6 


2.1/2.0 


450 < Q 2 < 1000 


2.04 


2.9 


4.2 


2.2/0.4 


1.94 


3.0 


5.4 


2.4/2.5 


1000 < Q 2 < 2000 


1.97 


5.5 


9.0 


2.5/0.2 


1.65 


6.2 


6.9 


2.5/2.7 


2000 < Q 2 < 8000 


1.92 


9.6 


11.1 


1.4/0.3 


1.50 


11.1 


14.4 


4.6/5.1 




0.3 < x 


p < 0.4 


100 < Q 2 < 175 


0.910 


4.0 


5.0 


1.1/1.9 


0.832 


4.1 


6.1 


1.5/2.2 


175 < Q 2 < 250 


0.925 


2.7 


4.0 


2.8/3.7 


0.892 


2.7 


4.0 


3.1/2.6 


250 < Q 2 < 450 


0.955 


2.6 


4.5 


3.9/2.3 


0.856 


2.8 


4.0 


3.0/2.9 


450 < Q 2 < 1000 


0.943 


4.2 


5.0 


2.9/3.0 


0.713 


4.9 


5.9 


1.8/1.8 


1000 < Q 2 < 2000 


0.807 


8.8 


11.3 


0.3/4.6 


0.691 


9.6 


11.1 


3.7/3.5 


2000 < Q 2 < 8000 


0.816 


14.2 


22.5 


1.9/2.6 


0.551 


18.2 


22.0 


1.5/3.4 




0.4 < x p < 0.5 


100 < Q 2 < 175 


0.453 


5.7 


8.3 


2.9/1.8 


0.395 


5.9 


7.1 


2.0/2.0 


175 < Q 2 < 250 


0.471 


3.7 


5.6 


4.9/3.1 


0.408 


4.0 


6.7 


3.1/3.8 


250 < Q 2 < 450 


0.437 


3.9 


6.0 


4.2/5.1 


0.390 


4.1 


5.2 


4.8/3.6 


450 < Q 2 < 1000 


0.505 


5.9 


9.0 


0.8/3.0 


0.346 


7.0 


9.1 


4.8/5.0 


1000 < Q 2 < 2000 


0.402 


12.6 


13.4 


2.4/1.1 


0.253 


15.3 


19.2 


4.9/3.8 


2000 < Q 2 < 8000 


0.445 


19.9 


21.3 


4.0/1.5 


0.216 


30.4 


37.5 


5.0/1.4 




0.5 < x 


p < 0.7 


100 < Q 2 < 175 


0.171 


6.4 


11.1 


3.3/2.3 


0.164 


6.7 


7.9 


3.5/1.3 


175 < Q 2 < 250 


0.191 


4.2 


6.3 


3.5/6.3 


0.146 


4.6 


5.9 


6.6/6.9 


250 < Q 2 < 450 


0.178 


4.4 


6.6 


4.8/4.9 


0.148 


4.8 


6.1 


6.7/5.6 


450 < Q 2 < 1000 


0.181 


6.7 


7.8 


4.3/6.2 


0.138 


7.9 


11.0 


7.0/3.3 


1000 < Q 2 < 2000 


0.143 


14.2 


15.9 


4.3/4.7 


0.106 


16.6 


17.5 


3.7/1.9 


2000 < Q 2 < 8000 


0.194 


21.7 


24.2 


9.0/3.6 


0.104 


29.1 


29.5 


0.8/5.5 




0.7 < x 


p < 1.0 


100 < Q 2 < 175 


0.0324 


12.2 


13.1 


3.8/5.1 


0.0297 


12.9 


15.8 


4.0/5.8 


175 < Q 2 < 250 


0.0332 


7.5 


8.6 


12.1/10.0 


0.0240 


9.0 


10.2 


13.4/14.5 


250 < Q 2 < 450 


0.0337 


8.2 


8.9 


9.6/7.5 


0.0216 


9.4 


10.4 


8.6/8.4 


450 < Q 2 < 1000 


0.0329 


12.9 


16.9 


7.2/1.8 


0.0200 


15.8 


19.9 


4.6/12.7 


1000 < Q 2 < 2000 


0.0274 


26.4 


51.0 


1.0/5.0 


0.0102 


39.8 


42.6 


5.1/11.5 


2000 < Q 2 < 8000 


0.0181 


43.6 


48.2 


13.7/8.3 


0.0083 


60.9 


64.4 


3.4/6.6 



Table 4: The measured normalised distribution of the scaled momentum for positively, 
D + (x p , Q), and negatively, D~(x p , Q), charged particles, as a function of Q 2 for different x p 
intervals, shown with the statistical error (<5 stat ), the total error including statistical and system- 
atic errors added in quadrature (5 to t ), and the correlated error coming from the electron energy 
scale uncertainty (5 sca [ e ) which is shown as two numbers (+/— ) and is not included in the total 
error. 13 



Q 2 [GcV^] 


A(i p ,Q) | 6 s tat *tot 




0.0 < x p < 0.02 


inn • f~\^ 1 

IUU < V < 1/0 


0.007 


0.036 


0.045 


± f D <. <. DU 


0.005 


0.019 


0.032 


iiOU <^ <. 4oU 


0.010 


0.015 


0.030 


4oU <. <. IUUU 


0.021 


0.016 


0.032 


1000 < Q z < 2000 


-0.019 


0.022 


0.034 


2000 < Q' 2 < 8000 


0.028 


0.027 


0.039 




0.02 < x p < 0.05 


IUU <. (tj <. 1 / J 


-0.008 


0.015 


0.030 


1 7e O-^ OKfl 

± / D <. Is^ <. ^JOU 


0.013 


0.009 


0.027 


iiDU <. <. 4oU 


0.003 


0.008 


0.028 


4oU <. <. IUUU 


0.012 


0.010 


0.028 


1000 < < 2000 


0.001 


0.017 


0.033 


2000 < Q 2 < 8000 


0.017 


0.027 


0.039 




0.05 < x p < 0.1 


IUU <. (tj < 1 / J 


0.010 


0.012 


0.028 


ll£> < < ZOU 


0.008 


0.008 


0.028 


ZOV <i 4bU 


0.003 


0.007 


0.026 


4oU <. <. IUUU 


0.011 


0.011 


0.029 


1000 < Q 2 < 2000 


0.003 


0.020 


0.033 


2000 < Q 2 < 8000 


0.004 


0.034 


0.046 




0.1 < Xp < 0.2 


iuu < y < 1/0 


0.008 


0.012 


0.029 


1|0 < <. zou 


0.001 


0.008 


0.027 


iiOU <. (4^ <. 4DU 


0.020 


0.008 


0.027 


4DU <L ^ IUUU 


0.003 


0.012 


0.028 


1000 < < 2000 


0.025 


0.024 


0.042 


2000 < Q 2 < 8000 


0.056 


0.041 


0.051 




0.2 < x p < 0.3 


1 nn ^ ^ 1 7r 
IUU <. (tj < 1 / J 


0.020 


0.019 


0.032 


17c ^- f\£ ^- e n 


0.020 


0.013 


0.030 


-iOU <. 4DU 


0.028 


0.013 


0.029 


4oU <. v ^ IUUU 


0.043 


0.021 


0.041 


1000 < Q 2 < 2000 


0.12 


0.041 


0.059 


2000 < Q 2 < 8000 


0.18 


0.074 


0.11 




0.3 < x p < 0.4 


IUU < < 1/0 


0.058 


0.029 


0.041 


17c ^- oen 

IfO s. y <. DU 


0.018 


0.019 


0.032 


iiOU <. (4^ <. 4DU 


0.049 


0.019 


0.034 


4oU <. <. IUUU 


0.15 


0.032 


0.045 


1000 < Q 2 < 2000 


0.047 


0.062 


0.083 


2000 < Q 2 < 8000 


0.32 


0.12 


0.18 




0.4 < xp < 0.5 


1 nn ^ /"i— ^ 1 
IUU <. <. 1 / 


0.033 


0.041 


0.001 


IfO < <. ZOU 


0.11 


0.028 


0.056 


zOU < C,/ < 4oU 


0.036 


0.028 


0.045 


4oU <. <. IUUU 


0.13 


0.043 


0.082 


1000 < < 2000 


0.18 


0.094 


0.12 


2000 < Q 2 < 8000 


0.39 


0.15 


0.18 




0.5 < xp < 0.7 


inn /"I 17c 
IUU < < IfO 


0.070 


0.047 


0.074 




0.12 


0.031 


0.048 


ncn ^- Hen 


0.081 


0.032 


0.047 


ARn d- 2 - <r 1 nnn 


0.095 


0.048 


0.068 


1000 < Q' z < 2000 


0.19 


0.11 


0.12 


2000 < Q z < 8000 


0.30 


0.17 


0.19 




0.7 < x p < 1.0 


100 < Q- 1 < 175 


0.014 


0.084 


0.096 


175 < Q' 2 < 250 


0.17 


0.055 


0.007 


250 < Q' 2 < 450 


0.23 


0.059 


0.068 


450 < Q 2 - < 1000 


0.25 


0.10 


0.12 


1000 < Q 2 < 2000 


0.33 


0.19 


0.25 


2000 < Q' 2 < 8000 


0.62 


0.40 


0.52 



Table 5: The charge asymmetry, A(x p , Q), as a function of Q 2 for different x p intervals shown 
with the statistical error (5 sta t) and the total error including statistical and systematic errors 
added in quadrature (8 to t)- 
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Figure 1: (a) The measured normalised distributions of the scaled momentum, D(x p ), for all 
charged particles and for positively (pos), D + (x p ), and negatively (neg), D~(x p ), charged par- 
ticles, and (b, c) the charge asymmetry, A(x p ), as a function of x p . The error bars include 
statistical (inner), and statistical plus systematic errors added in quadrature (outer). The data 
are compared to predictions from different models of the parton cascade and hadronisation pro- 
cesses implemented in leading order matrix element Monte Carlo programs and to the parton 
level before hadronisation. 
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Figure 2: The measured normalised distributions of the scaled momentum, D(x p , Q), for all 
charged particles and for positively (pos), D + (x p , Q), and negatively (neg), D~(x p , Q), charged 
particles separately, as a function of Q for nine different x p regions. The error bars include 
statistical (inner), and statistical plus systematic errors added in quadrature (outer). The data 
are displayed at the average value of Q. The data are compared to predictions from the CDM 
Monte Carlo program. 
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Figure 3: The charge asymmetry, A(x p , Q), as a function of Q for nine different x p regions. The 
error bars include statistical (inner), and statistical plus systematic errors added in quadrature 
(outer). The data are displayed at the average value of Q. The data are compared to predictions 
from different models of the parton cascade and hadronisation processes implemented in leading 
order matrix element Monte Carlo programs. 
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